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ABSTRACT 


Modifications  were  made  to  a  motion  prediction 
computer  program  for  a  hull  borne  hydrofoil  in  waves 
to  include  the  zero  speed  case.  The  program  can  now 
be  used  to  compute  the  si x-degrees-of- freedom  (6D0F) 
hull  borne  hydrofoil  craft  motions  for  a  craft  speed 
of  zero  as  well  as  for  a  constant  forward  speed,  less 
than  the  critical  "lift-off"  speed,  with  arbitrary 
heading  in  regular  waves.  Basically,  the  program 
consists  of  the  already  existing  "DTNSRDC  (David  W. 

Taylor  Naval  Ship  Research  and  Development  Center) 

6D0F  Ship-Motion  and  Sea-Load  Computer  Program"!*, 
which  computes  the  motions  of  hulls  in  the  dis¬ 
placement  mode,  with  appropriate  modifications  to 
incorporate  the  effects  of  the  submerged  foil  system 
on  craft  motions.  Presented  in  this  report  are  a 
discussion  of  the  mathematical  model  for  the  zero 
craft  speed  case,  input  information  pertaining  to 
the  hydrofoil  system,  and  a  discussion  of  the  results. 

The  Appendices  present  a  listing  of  the  program 
update  cords  and  the  subroutines  pertaining  to 
the  foil  system. 

ADMINISTRATIVE  INFORMATION 

This  work  has  been  authorized  by  and  sponsored  by  the  Naval  Material 
Command  (08T),  funded  under  the  Ship,  Subs,  and  Boats  Program,  Element 
62543N,  Task  Area  ZF-43-421-001  and  administered  by  the  Ship  Performance 
Department  High  Performance  Vehicles  Program  (1507)  under  Work  Unit  Number 
1500-100. 


INTRODUCTION 

A  description  of  the  "DTNSRDC  Hullborne  Hydrofoil  6D0F  Motion 
Prediction  Computer  Program"  for  the  zero  craft  speed  case  is  presented 
in  this  report.  It  serves  as  an  extension  to  the  previously  developed 
program  in  Reference  2,  which  predicts  the  motion  for  a  hullborne 
hydrofoil  craft  in  six-degrees-of- freedom  (6D0F)  advancing  at  a  constant 
forward  speed  in  the  displacement  mode  with  submerged  foils  at  an 

♦References  listed  on  Page  31. 
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arbitrary  heading  in  regular  waves.  The  program  is  an  adaption  of  the 
already  existing  "DTNSRDC  6D0F  Ship-Motion  and  Sea-Load  Computer  Program"*, 
based  on  the  theory  by  Salvesen,  Tuck,  and  Faltinsen'*,  which  was  developed 
for  the  prediction  of  motions  and  dynamic  loads  of  conventional  displace¬ 
ment  type  hulls,  and  is  utilized  for  planing  hulls  in  the  displacement 
mode  as  well.  The  program  modifications  consist  of  adding  the  hydrodynam¬ 
ic  force  terms  attributed  to  the  immersed  hydrofoil  system  to  the  already 
existing  equations  of  motion  for  the  hull.  This  technique  has  been  used 
successfully  by  R.T.Schmitke  who  developed  two  computer  programs  for 
simulating  the  motions  of  a  hullborne  hydrofoil  vessel  at  zero  speed  as 

well  as  for  forward  speeds  less  than  the  critical  "lift-off"  speed  for 

4  5 

head  seas  and  beam  seas  . 

The  hull  related  input  information  for  the  herein  presented  program 
remains  identical  to  the  original,  "DTNSRDC  6D0F  Ship-Motion  and  Sea-Load 
Computer  Program."  The  additional  input  information  relating  the  hydro¬ 
foil  system  is  given  for  computing  the  hullborne  hydrofoil  motions  either 
at  zero  speed  or  at  a  forward  speed  less  than  the  critical  "lift-off" 
speed  for  an  arbitrary  heading. 

The  modified  program's  output  consists  of  the  amplitudes  and  phases 
in  surge,  sway,  heave,  roll,  pitch  and  yaw  for  a  given  set  of  wave  fre¬ 
quencies  and  a  specified  set  of  forward  ship  speeds  and  headings,  as  well 
as  the  corresponding  response  amplitude  operators,  nondimensional  transfer 
functions,  and  nondimensionalized  excitation  forces  and  moments  as  in  the 
original  program  for  ship  motions.  Also  printed  are  the  coefficients  of 
motion  and  the  excitation  forces  for  the  hydrofoil  system  as  well  as  for 
the  hull.  Optionally,  one  can  obtain  the  two  sets  of  coupled  differential 


2 


equations  of  motion  in  matrix  form;  one  for  surge,  heave,  and  pitch  and 
the  other  for  sway,  roll,  and  yaw.  Each  set  is  given  for  the  hull  portion, 
the  foil-strut  portion,  the  combination  of  the  two,  and  the  final  inverted 
matrix  with  the  solutions.  Both  sets  are  for  the  minimum  specified 
frequency. 

MATHEMATICAL  MODEL 

The  analytical  model  for  determining  the  motions  of  a  hull  borne  hydro¬ 
foil  craft  at  zero  speed  is  obtained  by  adding  linear  and  quasi-linear 
terms  to  the  hull  terms  obtained  by  strip  theory  which  is  based  on  the 
solution  of  the  sectional  two-dimensional  problem  utilizing  a  close-fit 
source-distribution  technique  developed  by  W.  Frank^.  The  major 
assumptions  and  limitations  are: 

1)  wave  amplitude  of  unidirectional  regular  waves  and  craft  dis¬ 
placements  from  equilibrium  are  both  small. 

2)  both  the  craft's  beam  and  draft  are  much  smaller  than  its  length 

3)  the  craft,  both  hull  and  hydrofoil  system,  are  laterally 
symmetrical 

4)  dynamic  lift  attributable  to  the  hull's  planing  surfaces  is 
insignificant 

5)  the  foil  system  is  divisible  and  representable  by  a  set  of 
rectangular  foil  elements 

6)  interaction  between  the  hull  and  hydrofoil  system  is  negligible 
as  is  the  interaction  between  the  foil  elements 

7)  the  hydrofoil  system's  contribution  to  surging  is  negligible. 

HULL  EQUATIONS  OF  MOTION 

The  conventions  used  in  the  hullborne  hydrofoil  craft  motion  program 
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are  the  same  as  the  "DTNSRDC  Ship-Motion  and  Sea-Load  Computer  Program." 

3 

The  following  will  briefly  restate  the  definitions  used  .  As  shown  in 
Figure  1,  the  vessel-oriented,  right-handed  coordinate  system  is  defined 
to  have  its  origin  in  the  plane  of  the  undisturbed  free  water  surface. 

The  positive  x  axis  passes  through  the  craft's  stern,  the  y  axis  is  to 
starboard,  and  the  z  axis  is  vertically  upward  passing  through  the  center 
of  gravity.  The  heading  angle  y  is  defined  to  be  0  degrees  for  waves 
approaching  the  stern  and  180  degrees  for  waves  approaching  the  bow  as 
illustrated  in  Figure  2.  The  encounter  frequency,  we>  being  equal  to  the 
circular  wave  frequency  w  for  zer  craft  speed  is  given  by  the  expression 

a)  =  /2irg  ,  . 

X  (1) 

for  depths  much  greater  than  the  wavelength  A. 

With  the  assumption  that  the  equations  of  motion  are  either  linear 

or  have  been  linearized,  motion  displacements  can  be  expressed  as 

^  ^ j cos  ( cut—  g j ) ,  j  —  1 , . .  . , 6  (2) 

where  a.  are  the  amplitudes  and  e.  are  the  phase  lags  of  the  motion  with 
J  J 

respect  to  the  maximum  wave  elevation  above  the  origin.  The  subscripts 
j  =  1...6,  refer  respectively  to  the  translatory  displacements  of  surge, 
sway,  and  heave  and  the  angular  displacements  of  roll,  pitch,  and  yaw. 
Following  from  the  above  assumptions,  the  six  linear  coupled  differential 
equations  of  motion  can  be  written  in  complex  form  as: 

j.  (”jk  *  V\  *  BJk\  +  C,knk  -  F  e'V;  j  .  1 . 6  (3) 


where  m^  are  the  components  of  the  craft's  generalized  mass  matrix, 

are  the  added-mass  coefficients,  B.,  and  C..  are  the  complex  damping  and 

J  K  J  K 
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restoring  coefficients  and  F.  are  the  complex  amplitudes  of  the  excitation 
forces  and  moments. 

For  a  hullborne  hydrofoil  craft  as  well  as  for  a  conventional  dis¬ 
placement  hull  the  six  coupled  equations  of  motions  can  be  separated  into 
two  sets  of  equations.  With  the  exclusion  of  hydrostatic  restoring  co¬ 
efficients  that  are  equal  to  zero  for  both  the  hull  and  foil  system,  the 
first  set  of  three  coupled  equations  of  motion  in  surge,  heave,  and  pitch 
is 


Surge  (An  +  m)„j  +  B^  +  +  A15n5  + 


=  Fle' 


(4) 


Heave  A^r^  +  B^rij  +  (A33  +  m)n3  +  B33n3  +  C33n3  +  A35n5 


+  B35n5  +  C35n5  F3e 


(5) 


Pitch  A^nj  +  Bgjnj  +  A53n3  +  B53n3  +  C53n3  +  (A55  +  I5)n5 

'5 


+  B55^5  +  C55n5  *  F&e  e 


and  the  second  set  of  equations  is 


Sway 


Roll 


(^22  m^n2  B22n2  ^24  ”  ^G^4  B24^4 

+  A26n6  +  B26n6  +  C26n6  =  F2e 

(^42  “  mZg)n2  +  ^42^2  +  ^44  +  ^4^4  +  ^B44 
+^44n4  +  (A46  '  *46^6  +  B46n6  +  C46n6  =  F4e 


+  C24n4 

-itOgt 


+  ^44^4 

"V 


(7) 

(8) 
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Yaw 


••  ^  •• 

A62n2  +  B6292  +  ^A64  "  Wn4  +  B64^4  +  C64n4 


-i(o  t 


+  (A66  +  I6)n6  +  B66^6  +  C66n6  =  F6e’  “e 


(9) 


where  m  is  the  vessel's  mass,  Zg  is  the  location  of  the  vertical  center 
of  gravity  on  the  z  axis,  I.  is  the  moment  of  inertia  in  the  jth  mode, 

J 

and  I.,  is  the  product  of  inertia. 

J  K 

Considering  the  hull  alone,  as  in  the  hydrofoil  retracted  mode,  B*^ 
in  the  roll  equation  is  the  nonlinear  viscous  damping  of  the  hull.  It  is 
in  the  form  of  a  quasi-linear  function  in  terms  of  w,  viscosity,  hull 
geometry,  and  the  maximum  roll  amplitude  for  a  given  wave  slope  and 
heading,  at  zero  speed.  The  second  set  of  equations  in  sway,  roll,  and 
yaw  may  have  to  be  solved  a  number  of  times  until  the  difference  between 
the  maximum  estimated  roll  angle  and  the  computed  maximum  roll  angle  is 
within  an  acceptable  tolerance.  In  the  program  the  allowed  tolerance  is 
one  degree. 

The  derivation  of  the  hull  added  mass,  damping  and  hydrostatic  re¬ 
storing  coefficients  and  the  exciting  forces  and  moments  used  in  the 
"DTNSRDC  Ship-Motion  and  Sea-Load  Computer  Program"  are  presented  as 
three-dimensional  problem  of  ship-sections  oscillating  in  the  free  surface 
which  are  determined  by  a  close-fit  source-distribution  method  presented 
in  Reference  6. 

With  the  assumed  insignificance  of  the  hull  and  foil  system  inter¬ 
action,  the  hull  and  foil  system  contriubitton  to  the  added  mass,  damping 
and  restoring  coefficients  and  the  forcing  functions  are  simply  additive. 

For  example,  the  added  mass  coefficients  A.,  can  be  expressed  as 

J  K 

A  =  fiH  +  flF 

jk  Ajk  Ajk  (10) 
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H  F 

where  are  the  hull  added  masses  and  A^  are  the  frequency  independent 
foil  system  added  masses. 

HYDROFOIL  EQUATIONS  OF  MOTION 

In  the  derivation  of  the  hydrofoil  equations  of  motion  three  coor¬ 
dinate  systems  will  be  used.  One  system  is  the  inertial  reference  frame 
with  axes  xq,  yQ,  and  zq  and  unit  vectors  iQ,  j  ,  and  kQ,  respectively. 

The  second  coordinate  system  is  fixed  with  respect  to  the  hull  borne 
hydrofoil  craft  such  that  the  z  =  0  plane  lies  in  the  still  waterplane; 
z  axis  passes  through  the  craft's  CG;  the  positive  x  axis  passes  through 
the  stern;  and  the  positive  y  axis  points  to  starboard.  The  unit  vectors 
in  the  x,  y,  z  axes  are  i,  j,  k,  respectively.  The  third  coordinate 
system,  denoted  by  prime  ('),  is  also  fixed  with  respect  to  the  body. 

In  this  system  z '  =  0  is  the  equation  for  the  foil  element  with  the  origin 
(x  =  X,  y  =  Y,  z  =  Z)  at  its  midpoint,  with  the  axes  x  and  x'  parallel  to 
each  other.  The  primed  system  is  thus  translated  with  respect  to  the 
unprimed  system  by  X,  Y,  Z  and  rotated  about  the  x'  axis  by  the  dihedral 
angle  r,  as  shown  in  Figure  3. 

Hydrofoil  Coefficients 

The  force  on  a  rectangular  foil  element  with  chord  c  and  span  s  at 
zero  forward  speed  is  given  by  the  expression 

a(*N  +  V  +  Bi^vN  +  V2  +  b2VN  (1U) 

where  a($N  +  + 

a  =  added  mass  coefficient 
bj  =  viscous  damping  coefficient 
b^  =  wave-making  damping  coefficient 

v^  =  normal  velocity  component  on  the  foil  element  due  to  craft  motion 
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w^  =  normal  wave  velocity  component  on  the  foil  element 
Upon  linearizing  the  viscous  damping  term  by  the  principle  of  equivalent 
energy  balance  which  equates  the  energy  dissipated  by  a  non  linear  damped 


system  to  that  of  a  linear  damped  system,  equation  (11a)  becomes 
a(*N  +  V  +  3T  bll VN  +  Wn'(vN  +  V  +  b2vN 


(lib) 


In  determining  the  normal  velocity  v^  and  the  normal  acceleration 
v^  one  first  notes  that  the  velocity  of  the  foil  element  with  respect 
to  the  inertial  system  in  vector  notation  is 


v  =  v  +  (w  x  r) 

0  i  •  it  i  ; 

1  J  k 


.1 

n2 


04  05  Oe 

x  y  z 


With  the  unit  normal  on  t 


(12a) 

(12b) 


he  element  as 

n  =  -j  sin  r  +  k  cos  T  (13) 

the  normal  velocity  on  the  foil  element  in  the  body  reference  system 
becomes 

n)n  (14a) 


VN  =  <V0 


=  -sinr(n9  +  xnfi  -  7nd)  +  cosr(ri-  +  yfi.  +  xnK)  n 


(14b) 


and  the  normal  acceleration 

vN  =  -sinr(n2  +  xn5  -  zn4)  +  cosr(n3  +  yo4  -  xn5)  n  (15) 

For  harmonic  oscillations  with  a  frequency  of  u>,  the  normal  velocity  and 

acceleration  can  also  be  expressed  as 

VN  =  vN  elut  and  ^  =  i^  e1wt 

where  denotes  complex  amplitude.  Likewise, 

-  i  wt 

Oj  H  Oj  e  ;  J  =  2,  6 
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and  ni  =  i nift  • 

where  i  =  /T.  It  is  implied  that  only  the  real  part  on  the  right-hand 
sides  is  retained,  for  the  real  time  domain  representation. 

To  obtain  the  normal  wave  velocity  w^  and  acceleration  with  a 
wave  frequency  of  <y,  one  first  notes  the  complex  orbital  wave  velocity 
potential  at  a  point  x  ,  y  ,  z  in  the  inertial  reference  system  to  be 
<P  *=  -ic£  e  zQ  -  i  (yQ  siny  +  xQ  cosy)  +  iut  (16) 

2 

where  s  is  the  amplitude,  k  is  the  wave  number  u  / g,  and  y  is  the  heading 
angle  as  defined  and  illustrated  in  Figure  2.  The  normal  wave  velocity 
on  the  foil  element  at  its  mean  position  in  the  body  reference  system  is 
w  =  (v*  •  n)n 

N  =  Cw  -siny  sinr  +  i  cosr  exp[  k  z  -  i (y  siny 

+  x  cosy)  +  iut}i  n  (17) 

and  the  acceleration 

WN  =  luWN  (18) 

The  expressions  for  v^,  w^,  and  can  now  be  transformed  into 

the  primed  system  where 
x  =  X  +  x' 
y  =  Y  +  y'  cosr 
z  =  Z  +  y*  sinr 

(at  z’  =  0) 

With  the  neglect  of  chord-wise  variations  the  normal  velocities  in  the 
primed  system  become 

VN  "  (VN  +  y'  V"  09a) 

*  “  (V„  +  V*  n*)n  (19b) 
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where 


(19c) 

(I9d) 


WN  *  W^exp[iut  +  ky'(slnr  -  I  cosTsln  p)]  n 
wN  -  law' 

\  *  “Slnr(n2  +  Xr^  -  Zri^)  + 
cosr(n3  +  Yn^  -  x^) 

\  *  -sinr(n2  +  Xn6  -  Zn^)  + 
cosr(H3  +  Yn4  -  XH5) 

WN  =  Cw  (-sin  p  s i nr  +  i  cosr) 

• exp{k [Z  -  f  (Y  sin  p  +  X  cos  p)]} 


Inertial  Forces 

The  inertial  forces  on  the  foil  element  are  now  obtained  by  evaluating 

the  first  term  in  (lib);  using  the  accelerations  v^  and  w^,  using 
2 

0.25ttpc  for  a,  the  added  hydrodynamic  mass  for  a  plate  of  unit  length 
(Reference  7);  and  integrating  over  the  span.  The  resulting  translational 
modes  in  surge,  sway,  and  heave  are 


Surge 

A1  =  0 

(21a) 

Sway 

A2  «  Oj  (D^ s i n^r-D2s i nTcosT)  +02510? 

(21b) 

Heave 

A^  =  0j (DjSinrcosr-D2cos2r)  -  Q^osr 

(21c) 

where 

Qj  =  0.25irpc2s 

D,  ■=  n2  +  xH6  -  Z'r^ 

°2  “  S  +  "  X^5 

Qj  -  -0.257rpc2i>2  [cost  +  i  sin  p  slnT]  X 

s 1 nh [0. 5ks (s InT  -  i  sin  p  cosT)l 
k(sinr  -  7  s  In  p  cosT)  X 

•exp  { k [Z- I  (Ysinp  +  X  cos  p)] 

+  iut) 
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The  rotational  modes  in  roll,  pitch,  and  yaw  can  be  expressed  as  follows 


Roll 

\ 

“V 

-  V  + 

(2  Id) 

Pitch 

A5 

-  AjX 

+  A5 

(21e) 

Yaw 

A6 

X 

ts 

< 

1 

+  A6 

(21  f) 

where  the  last  terms  are  the  hydrodynamic  moments  about  the  applicable 


axis  passing  throuqh  the  element's  midpoint.  They  are 

Ai  ■  fl  *4 
.2 


Ai  -  o,  52  cos2  rn5 


5 

Ae  ’  “i  &S'"H 

Viscous  Damping  Forces 


The  linearized  viscous  damping  term  in  (lib)  with  b,  evaluated  for  a 
foil  element  of  unit  length  is 


£  pCDC  >VN  +  WN|(VN  +  V  \ 


(22) 


CQ  is  the  drag  coefficient  for  a  plate  with  an  angle  of  attack,  ct,  to  the 


direction  of  motion.  From  Reference  8 
0.0467a  a  <  40® 


1.17 


a  >  40° 


The  angle  of  attack  a'  is  given  for  0°  _<_  r  <  180°  as 


pitch  and  heave 
sway  and  yaw 
rol  1 


a 

a' 


1 90°-r I 


Y 

|arctan(  -■=■)-  T 


surge 

wherefrom 

a  ■  a' 

a  -  180° 


With 


0#  <  a'  1  90° 
90*  <_  a'  £  180’ 


»-  *  w-  -  (v-  +  »')  , 


I  tot 


(23) 
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p 


The  complex  amplitude  can  be  determined  as 

5n  +  ”n  ■  SN  +  1  “V' 


where 


+  exp[ky'(slnr  -  i  cosTsIny)] 


VN  "  1  t»[sinr  (n2  +  xn6  -  Zn^) 
+  cosr(n3  +  Yn4  -  x^)] 


(24) 


The  absolute  value  thereof  is 

lVN  +  WnI  +  [(vN  +  “n5  (vN  +  WN^*]1/2  (25) 

where  *  denotes  the  complex  conjugate.  Expanding  (25)  in  a  Taylor  Series 
(assuming  it  is  analytic  in  y1)  about  the  axis  y1  =  0  and  retaining  only 
the  first  order  terms  result  in 

(26) 


|vj|  +  w'|  -  f(y')  =  Hj  +  H2y' 

where  h  =  [V  V*  +  2R  (5  0  *)  +  W  W  *1 
1  1  N  N  1  N  N  ’  YU  J 


1/2 


H2  “  "if"  +  [l*-V  WNk(sinr  -  i  cosrstny)]*} 


with  R  referring  to  the  real  part  of  a  complex  quantity.  The  viscous 
forces  can  now  be  determined  by  substituting  (26)  and  (22)  and  integrating 


over  the  span,  that  is 

37  p  CDc  /  f(y')  (v^  +  w>)  dy'  (27) 

In  terms  of  the  degrees  of  freedom,  the  forces  are 

Surge  =  0  (28a) 

Sway  B2  ■  •—  pCpC  sin  T  P^  (28b) 

Heave  B^  *  yjj-  pC^c  cos  TP^  (28c) 

Roll  Bl»  "  37  PCDc  P2  (28d) 

Pitch  Bj  ■  pC^c  x  cos  rpj  (28e) 

Yaw  B6  "  37  Pcdc  X  5,n  rPl  (28f) 
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(29) 


where  Pj  and  P2  are  the  following  integrals 
P,  -  /f(y')  (v^  +  w^)  dy' 

•  »  1  tot — 

-  VlO  +  nlt  ail  *  WNe  °20 

p,  -  /f(y')  (vj  +  wjj)  (z  slur  +  v  cosr  +  y')dy' 

-  (zslnr  +  Ycosrjp,  +  VNQn  +  n4Q12  +  WNe 
The  subscripted  Q's  are  evaluated  as  follows: 

Ql0  "  /f(y')dy'  -  HjS 

H  s3 

Qn  “  /f(y')y'dy  - 

H  s3 

Q12  -  /f(y')y,2dy'  =  -L- 


020 

02! 


/f(y')eGy'dy' 

2("r  '  ^  sfnh(^-)  +  ^  cosh(^) 


;f(y')y'eGy  dy1 


**H9  2Hi  H,s  g 

(“j"  +  ~~2~  +  ~2G~^  slnh(T^  +  Cf  ^H1 

G  G 


and  G  =  k(sinr  -  siny  cosr) 


2H.  rc 

COSh(-y) 


Wave-making  Damping 

A  foil  element  intersecting  or  near  the  free  surface  will  generate 
waves  when  oscillated  laterally.  Although  this  resitance  is  of  lesser 
importance  it  should  nonetheless  be  included  for  the  zero  speed  case. 
The  wave-making  damping  for  a  vertical  foil  element  oscillating  in  sway 


is 


IT  -  C. 

-  puCwcs  n2 


(31) 


2 

where  C  is  a  function  of  w  s/g  (Reference  5).  A  plot  of  C  obtained  by 
w  w 

the  Frank  close  fit  method  is  given  in  Figure  4.  From  (34)  the  following 
damping  forces  have  been  obtained  for  an  element  with  a  dihedral  angle  of 
r 
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Surge 

C1 

-  0 

(32a) 

Sway 

C2 

■  f  PU)Cwcs2sInr  n2 

(32b) 

Heave 

C3 

-  J  pa)Cwcs  COST 

(32c) 

Roll 

C4 

-  |-p<aCwcs2(Z  si nT  +  Y  COST) 

(32d) 

Pitch 

C5 

■  j-  poises2  X  cosTn^ 

(32e) 

Yaw 

C6 

“  j  poises2  X  slnTrig 

(32f) 

The  hydrofoil  related  forces  and  moments  developed  and  given  in 

equation  (21),  (28),  and  (32)  as  A.,  B.,  C.  for  the  six  degrees  of  freedom 

*1  J  J 

j  =  1,  6  are  now  multiplied  by  a  correction  factor  for  finite  spans. 

This  factor  is 

where  AR  is  the  foil  element's  aspect  ratio  s/c.  The  contributions  of  all 

the  individual  foil  elements  of  the  complete  hydrofoil  system  are  then 

summed  and  the  resulting  total  of  A.,  B.,  and  C.  separable  into  the  motion 

J  J  J 

and  wave  excitation  terms,  are  added  to  the  corresponding  hull  related 
terms  in  Equations  4  through  9  with  the  assumption  that  additional  hull- 
foil  interaction  terms  are  negligible. 

Further  simplification  was  required  to  make  the  program  more  tract¬ 
able  by  neglecting  some  foil  system  attributable  cross  coupling  terms; 
justified  on  the  basis  that  their  contributions  were  insignificant  as 
evident  from  computed  and  experimental  results.  Also  simplified  was  the 
quadratic  damping  force  by  treating  v^  and  w^  independently  whereby 
equation  (22)  becomes 

37  p  cDclvNivN  and  37  p  CDCIWnIWN 
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COMPUTER  PROGRAM 


Based  on  the  foregoing  mathematical  model,  the  program,  as  presented 
in  Reference  2,  was  expanded  to  include  the  hull  borne  hydrofoil  motions 
in  regular  waves  at  zero  speed  as  well  as  forward  speed  at  an  arbitrary 
heading.  The  basic  program  was  the  existing  DTNSRDC  Ship-Motion  Computer 
Program.  In  itself  the  basic  program  can  determine  the  hydrofoil  craft's 
motions  in  the  foil  up  mode  in  6D0F  in  regular  unidirectional  waves  of  any 
heading.  The  modification  consists  of  adding  the  foil  system's  coefficients 
of  motion  and  its  excitation  forces  and  moments  to  the  corresponding  terms 
for  the  hull.  As  a  consequence,  three  card  sets  listed  below  and  pertaining 
to  the  foil  system  of  the  hull  borne  hydrofoil  craft  are  added  onto  the 
existing  34  data  card  sets. 

A.  Input  Description 

Data  Card  Set  35  -  one  card  with  format  (314) 

IFOIL:  2  for  hydrofoil  craft  in  the  foils  down  mode.  All  other 
integer  values  are  for  retracted  foil  systems  where  only 
the  hull  is  subjected  to  hydrodynamic  forces. 

NZERO:  "umber  of  zero  speed  cases  for  the  hull  borne  hydrofoil  in 
the  foils  down  mode 

IPRINT:  option  of  printing  the  matrix  equations  of  motion.  With 
I  PRINT  =  0  printing  of  matrices  is  suppressed  and  for 
IPRINT  =  1  printing  of  matrices  takes  place. 
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Data  Card  Set  36  -  One  card  with  format  (15,  3F12.2) 

NF:  total  number  of  foil  elements  on  the  starboard  side  of  the 

hydrofoil  craft.  This  total  consists  of  the  elements  in  symmetry 
about  the  xz-plane  plus  the  elements  lying  in  the  xz-plane. 

FVOL:  is  the  displaced  volume  of  the  entire  foil  system  (including 
the  portions  of  the  struts  that  are  immersed).  The  unit  is 
W0RD**3,  (see  Data  Card  Set  4  of  Reference  1) 

FXCB:  the  foil  system's  longitudinal  center  of  buoyancy,  LCB,  with 

respect  to  the  entire  craft's  LCB,  i.e.  the  x  value  in  the  body 
coordinate  system  in  units  of  WORD. 

FZCB:  the  foil  system's  vertical  center  of  buoyancy,  VCB,  i.e.  the  z 
value  in  the  body  coordinate  system  in  units  of  WORD. 

Data  Card  Set  37  -  one  card  per  foil  element  with  format  (F3.0,  5F7.2,  F5.0, 
F10.7,  F5.0,  F5.0) 

CPL:  If  the  plane  foil  element  lies  in  the  center  pic.n?  i.e  the  xz-plane, 
CPL  =  1  and  in  all  other  cases 

CPL  =  2  where  symmetry  of  the  elements  about  the  center  Diane 
is  assumed. 

The  remaining  parameters  in  this  set  pertain  to  the  foil  element  or  portion 
thereof  on  the  starboard  side  of  the  centerplane  or  for  the  element  lying 
in  the  centerplane. 

SPAN:  the  length  of  the  rectangular  foil  element  taken  in  a  line 
parallel  to  the  yz-plane  in  units  of  WORD. 

CHORD:  the  width  of  the  rectangular  foil  element  taken  in  a  line 
parallel  to  the  xz-plane  in  units  of  WORD. 

X:  x  coordinate  of  the  foil  element's  midpoint  in  units  of  WORD 
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Y:  y  coordinate  of  the  foil  element's  midpoint  in  units  of  WORD 
Z:  z  coordinate  of  the  foil  element's  midpoint  in  units  of  WORD 
DGAMMA:  is  the  dihedral  angle,  i.e.  the  angle  between  the  foil  element 

and  the  y-axis.  If  the  foil  lies  in  the  centerplane  DGAMMA  =  90. 
The  range  of  DGAMMA  is  between  0  and  180  degrees. 

CLZ:  is  the  vertical  lift  slope  of  the  foil  element  in  dimensionless 
units  (not  used  in  zero  speed  case) 

ASP:  is  a  positive  number  utilized  in  the  aspect  ratio  correction 
factor  AR/(AR  +  ASP)  for  foil  elements  of  finite  span.  For 
the  example  cited  ASP  =  3. 

STRUT:  0.  for  a  non-surface  piercing  foil  element.  All  other  values 
are  for  surface  piercing  foils.  This  control  variable  for 
the  wave-generated  damping  is  used  in  the  zero  speed  case 
only. 

Data  Card  Set  38  -  NZERO  cards,  one  for  each  hull  borne  hydrofoil  zero  speed 
case  foils  down  mode.  The  format  is  (2F8.2,  2F8.4) 

AMPM0T(2):  Estimated  maximum  sway  amplitude  in  units  of  WORD 
AMPM0T(3):  Estimated  maximum  heave  amplitude  in  units  of  WORD 
AMPM0T(5):  Estimated  maximum  pitch  amplitude  in  radians 
AMPM0T(6):  Estimated  maximum  yaw  amplitude  in  radians 

The  program  for  the  zero  speed  case  uses  a  "trial  and  error"  procedure 
for  solving  the  quasi-linear  set  of  equations  in  sway,  roll,  and  yaw.  Keeping 
the  waveslope  constant,  the  initial  maximum  sway,  roll,  and  yaw  amplitudes 
are  used  in  computing  the  corresponding  motion  responses.  The  program  then 
compares  the  maximum  computed  amplitudes  with  the  estimated  values  and  if 
the  differences  in  sway  or  yaw  is  greater  than  8  percent  or  the  difference 
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in  roll  is  greater  than  1  degree  new  estimates  are  determined  where  needed 

and  used  in  the  repeated  computation  for  motion  responses.  These  iterations 
take  place  until  the  computed  and  estimated  maximum  amplitudes  in  sway, 

roll,  and  yaw  fall  within  the  acceptable  limits.  If  the  proper  closure 
still  does  not  occur  for  the  fifth  iteration,  the  program  prints  the  results 
for  the  last  iteration  which  includes  a  table  of  the  estimated  and  computed 
values  and  the  differences  and  a  note  as  to  the  response  where  the  proper 
closure  failed.  Although  it  is  advisable  to  have  good  estimates  in  maxi¬ 
mum  pitch  and  heave  amplitudes,  errors  in  these  are  not  as  critical.  As  a 
consequence  the  program  does  not  determine  new  estimates  for  pitch  arid  heave. 

The  program's  output  is  essentially  the  same  as  for  the  original  DTNSRDC 
Ship-Motion  Computer  Program  (Ref.l).  In  addition,  a  listing  of  the  foil 
related  input  cards,  Data  Card  Sets  35  through  38,  are  given  as  well  as  the 
foil  related  non-dimensi onal ized  coefficients  and  forcing  functions  and  a 
modified  listing  of  the  estimated  and  computed  maximum  response  amplitudes 
to  include  sway  and  yaw  as  well  as  roll.  As  an  option  the  matrix  equations 
for  the  hull  alone,  for  the  hydrofoil  system,  for  the  combined  hull  and  foil 
system,  and  the  solved  simultaneous  differential  equations  with  the  inverted 
matrix  can  be  printed. 

In  addition  to  the  original  subroutines  of  the  DTNSRDC  Ship-Motion  are 
the  subroutines  FOIL,  THEO,  and  EXCIT  for  the  hullborne  hydrofoil,  foils 
down  mode,  at  speeds  greater  than  zero  and  the  subroutine  ZERO  for  the  zero 
speed  case.  Appendix  A  presents  the  update  cards  needed  for  the  original 
DTNSRDC  Ship-Motion  Computer  Program  and  Appendix  B  presents  the  subroutines 
ZERO,  FOIL,  THEO,  and  EXCIT.  The  latter  two  subroutines  require  a  subprogram 
for  computing  Bessel  functions.  The  entire  program  is  on  digital  7-track 
magnetic  tape  which  can  be  used  on  the  CDC  6000  computer  system. 
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COMPARISON  OF  PREDICTED  AND  EXPERIMENTAL  RESULTS 

The  first  experimental  results  chosen  for  comparison  with  predicted 
results  were  the  hullborne  motion  measurements  of  the  313-ton  Plainview 
AG(EH)-1  hydrofoil  craft.  The  experiments  were  conducted  bv  R.M.  Vuolo  on 
a  1:12  scale  model  in  DTNSRDC's  Harold  E.  Saunders  Seakeeninn  and  Maneuverinn 
Facility.  The  model  was  run  in  both  the  foils  up  and  foils  down  modes  in 
unidirectional  regular  waves.  The  full  scale  speeds  were  6  and  12  knots  at 
the  three  headings  of  head  (180°),  bow  (150°),  and  beam  waves  (90°).  The 
regular  waves  were  of  constant  wave  steepness  1/60,  and  wave  lenqths  ranged 
from  L/X  =  0.33  to  L/X  =  4.0,  corresponding  to  wave  frequencies  of  w  =  0.57 
to  1.98  rad/sec. 

The  Hullborne  Hydrofoil  Six-Degree-of-Freedom  Motion  Program  was  like¬ 
wise  run  in  both  the  foils  up  and  foils  down  modes  at  the  three  headings  of 
180°,  150°  and  90°.  The  predicted  motions  at  craft  speeds  of  6  and  12  knots 
are  presented  and  compared  with  the  experimental  results  in  an  earlier 
report.  Reference  2.  At  zero  speed  the  predicted  motions  at  the  three 
headings  generally  agreed  well  with  the  experimental  results  as  nivne  in 
Figures  5a,  5b,  and  5c  which  show  the  craft's  transfer  functions  in  heave, 
roll,  and  pitch  versus  wave  frequency  together  with  the  phase  lag  with 
respect  to  the  maximum  height  of  the  wave  ut  the  CG.  The  effect  of  the 
immersed  hydrofoil  system  on  pitch  and  heave  is  insignificant  at  all  three 
headings  both  theoretically  and  experimentally.  Roll  is  the  only  deoree 
of  freedom  notably  affected  by  the  immersion  of  the  hydrofoil  system  as 
seen  in  the  90  and  150  degree  headings.  Figures  5b  and  5c.  Some  of  the 
minor  discrepancies  between  the  predicted  and  experimenal  results, 
especially  in  roll,  are  attributable  most  likely  to  inaccurate  roll 
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gyradii  estimates  since  the  model's  roll  gyradii  for  the  foil  up  and  the 
foil  down  modes  were  not  available.  A  less  significant  source  of  error 
may  also  be  inaccuracies  in  the  foil  system's  disDlacement  and  center  of 
buoyancy  estimates. 

The  second  set  of  experimental  results  used  for  comparison  with  the 
predicted  motions  of  the  Hull  borne  Hydrofoil  6D0F  Motion  Prediction  Pro¬ 
gram  was  a  1:20  scale  model  of  the  220  L-TON  PHM  hydrofoil  craft.  The 
mdoel  was  tested  at  the  Davidson  Laboratory  in  simulated  beam  Sea  States 
3  and  5  in  the  foils  submerged,  hull  borne  model  at  zero  speed,  Reference 
9.  Shown  in  Figure  6  are  the  measured  roll  responses  in  model  scale 
together  wi tn  predicted  roll  responses  which  agree  satisfactorily. 
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CONCLUDING  REMARKS 


The  motions  predicted  for  the  AG(EH)-1  hydrofoil  craft  at  zero  speed 
in  the  hull  borne  model  with  retracted  and  submerged  hydrofoils  are  in 
satisfactory  agreement  with  the  available  expeirmental  data  for  regular 
waves  at  90,  150  and  180  degrees  headings.  The  DTNSRDC  Hull  borne  Hydro¬ 
foil  bDOF  Motion  Prediction  Computer  Program  also  predicted  satisfactorily 
the  roll  response  at  zero  sneed  of  the  PHM  hydrofoil  craft  in  beam  Sea 
States  3  and  5  as  compared  to  experimental  results. 

Additional  comparisons  directed  toward  verifications  of  the  DTNSRDC 
Hull  borne  Hydrofoil  6D0F  Motion  Prediction  Computer  Proaram  should  be 
made  as  more  experimental  results  become  available.  Although  the  proaram 
can  be  used  at  present  to  predict  craft  motions  at  zero  and  forward  speeds, 
no  verification  has  been  made  of  craft  responses  at  wave  headings  less 
than  90  degrees  due  to  a  lack  of  experimental  data. 
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TRAN SLAT ORY  DISPLACEMENTS 


nj  -  SURGE 
U2  «  SWAY 
n  3  =  HEAVE 


ANGULAR  DISPLACEMENTS 


Ui+  =  ROLL 
Us  "  PITCH 
Ug  *  YAW 


Figure  1  -  Sign  Convention  of  Body  Coordinate  System 
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eading  Angle, 


Figure  3  -  Primed  Coordinate  System  and  Dihedral  Angle,  r, 
in  xyz  Body  Coordinate  System 
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(u)  circular  wave  frequency 
s  span 

g  gravitational  acceleration 


Figure  4  -  Strut  Wave-Making  Damping  Coefficient,  C 


Figure  5b  -  Non-Dimensional  Transfer  Function  and  Phase  Versus  ui 
of  the  AG(EH)- 1  in  Bow  Seas,  y=150° 
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Figure  5c  -  Non-Dimensional  Transfer  Function  and  Phase  Versus  to 
of  the  AG(EH)-1  in  Beam  Seas.  y=90° 
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Figure  6  -  -loll  Gain  versus  Wave  Frequency  of  the  PH*1  Hydrofoil 
Craft  In  Ream  Sea  States  3  and  5 
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APPENDIX  A 


LIST  OF  UPDATE  CARDS  FOR  DTNSRDC 
6D0F  SHIP-MOTION  PREDICTION  COMPUTER  PROGRAM 

The  update  cards  listed  below  are  required  to  modify  the  existing 

DTNSRDC  6D0F  Ship-Motion  and  Sea-Load  Computer  Program  to  the  Hullborne 

Hydrofoil  6D0F  Motion  Prediction  Computer  Program.  This  modification 

allows  the  use  of  the  hydrofoil  related  subroutines  presented  in  Appendix  B. 


*10  FMOO 

*  INSERT  LKO . 25 

COMMON  /P FOIL/  DUM5(165> 

COMMON  /QFOIL/  0UM6<36) 

*  INSERT  PR1.22 

COMMON  /PFOIL/  RHO, I F 01 L , N2ER0 , I PRI NT , NF , Q( 10,10) , AMPMOT C 6, 10 ) 

*0  PR1.255 

IF  C  I 0 AMP-2 )  9090,5055,9080 
•INSERT  PR1.28Q 

905  FORMAT)  I5.3F12.2) 

906  FORMAT(12X,I5,3F12.2) 

907  FORMAT (  I5,3F12.2J 

970  FORMAT {  F3. 0 , 5F7. 2 , F 5 . 0, FI 0 .7 , F5 . 1 , F5. 0) 

971  FORMAT(12X,F3.0,5F7.2,r5.0»F10.7*F5«i»F5»0) 

980  F0RMAT(//13H  CARD  SET  36) 

990  FORMAT (//13H  CARO  SET  37) 

7001  F ORMA T  <  2  14  ) 

7002  FORMAT ( / /22H  FOIL  OATA  CARO  INPUT) 

7003  FORMAT (/732HI  HYDROFOIL  VESSEL  KITH  FOILS  UP/) 

7004  FORMAT  t//34Hl  HYDROFOIL  VESSEL  KITH  FOILS  DOWN/) 

2212  FORMAT  (/21X,1H1, 9X, 1H2, 9X, 1H3,9X,1H4,9X, 1H5,9X, 1H6,9X,1H7,9X,1H8/ 
24X,8HC0LUMNS  , 8 ( 1 0H1 2 34 5 6 78 90 ) ) 

2213  FORMA T(/14X,2HNF,9X, 4 HF VOL , 8X , 4HF XCB, 8X,4HCZC8/) 

2214  FORMAT  </12X,3HCPL,3X,4HSPAN,2X, 54CH0R0, 3X, 4HX { S) , 4X, 1HY, 6X , 1H7, 2X 
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.....A. 


HL2CEDINQ  PAOX  BUW.KOT  IUMD 


ooooooooo  o  o  o  r>  o  o  o  anonnoo 


2,fcH0GAMMA,3X,3HCLZ,5X,3HASP/) 
DATA  CARO  SET  35 


HYDROFOIL  VESSEL  WITH  FOILS  UP 
HYOROFO I L  VESSEL  WITH  FOILS  DOWN 
PRINTOUT  OF  MATRIX  EQUATIONS  (NO 


-  IFOIL= 1 

-  IF  0IL  =  2 

=  0  ,  YES  =  i> 


RFADC5, 7001»  IFOIL, IPRINT 
IF ( IFOIL  .NE.  2)  IFOIL=l 
WRI TE (1 »  700 1)  IFOIL 
IF(IFOIL-l)  9091,9091,9092 

9091  WRITE (6, 7003) 

GO  TO  L515 

9092  WRITE<6, 70041 
WRITE (6,70021 
WRITE(6, 22121 


DATA  CARD  SET  36 


NUMBER  OF  INPUT  FOIL  ELEMENTS,  DISPLACED  VOLUME  (WORD**3), 
LONGITUDINAL  CENTER  OF  30YANCY  FROM  F. ?.  ANO  VERTICAL  CENTER  OF 
80YANCY  FROM  WATERLINE  (WORD)  OF  THE  ENTIRE  HYDROFOIL  SYSTEM 

REAO (5,905)  NF, F VOL,  FXCB , FZC8 
WRITE (6,980)  • 

WRITE(6, 22131 

WRI TE(6,906)  NF, F VOL , FXC 8 ,F ZCB 
WRITE(1,907)  NF,FV0L,FXC3,FZCB 
WRITE (6,9901 
WRITE (6,2214) 

DO  100  1=1, NF  ' 


DATA  CARD  SET  37 


FOIL  ELEMENT  IN  VERTICAL  CENTER  PLANE  (CPL=1.  FOR  YES,  CPL=2.  FOR 
NO)  ,  HYDROFOIL  ELEMENT  SPAN  (FT),  CHORD  (FT),  COORDINATES  X,Y,Z 
OF  MIOPOINT  (FT),  DIHEORAL  ANGLE  OF  V-F0IL  (DEG),  VERTICAL  LIFT 
SLOPE  ( W0RD2/W0RD ) ,  ASP  IN  THE  FACTOR  AR/(AR+ASP)  FOR  FINITE 
SPAN 

READ (5,970)  ( Q (N , I ) , N= 1 , 10 ) 

WRITE  (6,971)  (0(N,I)  , N= 1,  10) 

WRITE(1,970)  (Q(N,I) , N= 1 , 10) 

100  CONTINUE 

WPITE(6,2212) 

9515  CONTINUE 
*  INSERT  PR1.418 

IF(IFOIL-l)  51,51,52 
52  FXC0=FXC8/EL 
FZC8=FZCB/ELL 
FV0L=FV0L/EL3 
HVOL  =  TVOL 
TVOL=HVOL»FVOL 
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TPST=  (TPST*HVOL«-FXCB*FVOL)  /TVOL 
CBV  = (CBV*HVDL*F  ZCB*FVOL) / TVOL 
51  CONTINUE 

RHO^TMAS/ (TVOL* EL  31 
•I  SP5.21 

COMMON  /PFOIL/  RHO, I FOI L , NZERO , IPR INT , NF , Q (1 0 , 1 0 ) . AMPMOT 16, 10 ) 
COMMON  /QFOIL/  GA(6,6> 

*0  SP5.41 

DIMENSION  GMU(6,6),T0DE 1 5 , 6 ) , TEV F 1 6 , 6 > , BOOF ( 6 ) , BE VF (6 ) , T ( 32 , 62) 
DIMENSION  T00A(6,&!,TEVA(6,6),BDDA(6>,BEVA(b>,T0DB(6,6),TEVB(6»6)» 
2B0DBt  6)  ,BEVB(6)  ,  TOOC  (6,61  ,TEVC(S,6),300C(b),BEVC(6> 

*1  SP5.42 

BACKSPACE  1 

2198  FORMAT (118H1N0N-0IMENSI0NAL  ADOED  MASStDAM’INGjANO  RESTORING  COEFF 
2ICIENTS  ANO  EXCITING  FORCES  AND  MOMENTS  OF  ThE  STRUTS  AND  FOILS! 

2199  FORMAT (41X,12A&, 15X, 3H»**///17X, 9HHEADING  *,F5.0,4H  DEG  ,  7X , 1 2HSHIP 

2  SPEED  -,F6.2,6H  KNO T S/ 1 8 X , 16H (HE A D  SEAS  =1 80 ) , OX , 15HFR0UDE  NUMBER 

3  =,F7. !*,///) 

2200  F ORM AT {  /58H  NON-DIME NSIONALIZED  ADDED  HASS  COEFFICIENTS  OF  THE  FOI 
1LS-) 

22  01  F0RMAT(3X.6HHE(N0),5X»6HA(1,1)  ,6 X , 6H A ( 2 , 2 ) , bX  ,  5H A  (3,3),SX,6HA(4,4) 
1 , 6:.,6HA(5,5)  *6X,6HA(6,G>  ,6X,6HA<3,5>  ,  5  X ,  6  HA  ( 2 ,  S  )  ,  6  X,  6HA  (2,4),6X,6H 
2 A  (4,6)) 

2202  FORMAT (3X,F6.3,1P10E12.4) 

2203  FORMAT (9X.1P10E12. 4) 

2204  FORMAT  (  /  5  5  H  NON- D I  ME  N  SI  ON  AL I ZE  D  DAMPING  COEFFICIENTS  OF  THE  FOILS- 
1) 

2205  FORMA T(3X,6HWE(N0),5X,SHB(1 ,1) ,6X , 6HB ( 2 , 2 > , bX  ,  SHB (3,3) ,6X,6HB(4,4) 
1 , 6X.&HB (5,5) ,bX,6H3(6,6) , 6X ,6HB( 3,5) , 5X,&H3(2,3) ,bX,bHB  (2,4) ,6X,6H 
20  (4,6) ) 

2206  F0RMAT(/57H  NON-O IME NSI ON AL I  ZED  RESTORING  COEF-ICIENTS  OF  THE  FOIL 
1S-) 

2207  FORMAT (3X.6HWE (NO >,5X,6HC ( 1 , 1 ) , &  X , &HC I  2 , 2 > , 6X , 3 HC ( 3, 3 > , 6X , 6HC (4 , 4) 
1, 6X,6HC(5,5) ,6X,6HC(b,6) , 6X.6HCI 3,5) ,6X,6HG(2,5),6X,6HC(2,4),6X,6H 
2C  (4,6)) 

2208  FORMAT (/,*  NON-D I MENSION ALI ZED  FORCE  AND  MOMENT  FUNCTIONS  OF  THE  F 
10 ILS- * ) 

2209  FORMAT! 3X.6HWE (NO ) , 4 X ,5 HSURGE , 8X , 4HSN A Y , 7X, 5HHE A VE , 8X , 4HR0LL , 7X  ,  5H 
2PITCH.9X.3HYAW) 

2210  FORMAT (3X, F6. 3»1P6E12*4) 

2211  FORMAT (9X,1P6E12«4) 

8001  FORMAT (///51H  TOO*X=BOD  0  EF  ORE  INSERTION  OF  HYOROFOIL  ELEMENTS/) 

8011  FORMAT (///51H  TEV*X=BEV  BEFORE  INSERTION  OF  HYDROFOIL  ELEMENTS/) 

8003  FORMAT (///50H  T0D*X=80D  AFTER  INSERTION  OF  HYDROFOIL  ELEHENTS/J 

8013  FORMAT ( / //5 OH  TEV*X=BEV  AFTER  INSERTION  OF  HYDROFOIL  ELEMENTS/) 
6002  FORMAT!///*  MATRICI ES  TODF  AND  BO  DF  * / ) 

8012  FORMAT!///*  MATRICIES  TE VF  AND  BEV-*/> 

8021  FORMAT <9X, 1P6E12. 3,7X,1HX , 12, 7H  R E A L , 8 X , 1P1E 1 2. 3) 

8022  F0RMAT(9X,1P5E12.3,7X,1HX,I2,7H  IMA G , 8X , IPIE12. 3) 

8004  FORMAT ( ///14H  X  =  I NV T OD *  8 00/ > 

8014  FORMAT ( ///14H  X=  INV TEV * OE V/> 

8023  F0RMAT(9X,1PIE12.3,12X,1P6E12.3, 12X.1P1E12.3) 

*0  SP5.225 

IF ( IBILGE-1 )  3003,3003,3004 
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•  O  0*0  oooooo 


*0  SP5.266 

IF(IENO-l)  3001,3001,3002 
*1  SP5.306 

*****  **  **  •»» 

FOR  A  HULLB0RNE  HYDROFOIL  (IF0IL=2>,  SUBROUTINE  -FOIL-  CALCULATES 
THE  MOTION  COEFFICIENTS  ANO  THE  EXCITATION  FORCES  AND  MOMENTS  OUE 
TO  THE  FOILS.  RETURNEO  ARE  THE  TERMS  FOR  THE  TOO,  BOO,  TEV,  ANO 
BEV  M ATRICIES 

IF(IFOIL-l)  7200,7200,7100 
7100  WAMPL2=WAVAMP(LL) 

H0G2=H0G1 ( HM) 

7104  CALL  F0IHT0DF,TEVF,B0DF,8EVF,VKN0TS,RAMPL2,H0S2,GXI,ELL,RH0,NF,Q, 
2TMAS,T,LL,HVLNTHI 
7106  DO  7110  JA=1,6 
DO  7112  JB=1,6 
TOOA ( JA, J8)=T0D (JA, J81 
TEVA(JA,JB)=TEV(JA,JB) 

TOO ( JA, J8)=T0D( JA  ,J9)  *TODF( JA.JB) 

TEV (JA, J8»  =TEV< JA.JB) *TEVF< JA.JB I 
T009( JA.JB) =T00 (JA.JB) 

TEVB( JA.JB) =  TEV( JA.JB) 

7112  CONTINUE 
7110  CONTINUE 
7200  CONTINUE 

*****  **»»•*» 

I  SP5.418 

************ 

IF(IFOIL-l)  7600,7600,7500 
7500  DO  7510  J A= 1 , 6 

BOOA(JA)=BOO(JA,1) 

BEVA( JA)=BEV(JA,1) 

BOO(JA,l)=nOD(JA,l) *  BOOF ( JA I 
BEV ( J  A, 1)=9EV (JA, 1) ♦BEVFIJAJ 
B0DB(JA)=80D(JA,1) 

BEV8(JA)=BEV(JA,1> 

7510  CONTINUE 
7600  CONTINUE 

************ 

D  SPS.435 

DO  8201  IQ=1,6 
BOOC(IQ)=BOO(IQ,1> 

DO  8202  JQ=1,6 

8202  TODC(ia,JQ)=TOD(IQ,JQ) 

8201  CONTINUE 

IF(IO-l)  501,501,502 
* D  SP5.437 

DO  9203  IQ=1,6 
8EVC(IQ)=nEV(IQ,l) 

00  8204  JQ-l, 6 
8204  TEVC(IQ»JQ)=TEV(IQ,JQ) 

8203  CONTINUE  . 

IF(IO-l)  503,503,502 
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*0  SP5.461 

IF { IPRES-1 >  5202,5202,5203 
*  D  SP5.490 

IFIJS-l)  6222,6222,6223 
•D  SP5.531 

1505  IF ( IF  OIL- 1 )  2401,2401,2402 
2402  WRITE  16,2198) 

00  5614  JH= 1, 2 
IF ( JH  .EQ.  1)  H= 1 
I F ( JH  .EQ.  2)  H=  6 

IF  <  H  .EQ.  6  .AND.  PRNTOP  .EQ.  MIN)  GO  TO  5614 
WRITE (H, 2199)  TI TO, HD IG 1 , VKNOTS, F N( J J ) 

5614  CONTINUE 

WRITE(6,2200) 

WRITE (6,2201) 

00  2300  LX  =  1 , NOK 
LWEINC=NOK-LX*l 

2300  WRITE (6,2202)  ZN ( L WE  INC) , ( T (L WEI NC , KX) ,KX  =  1 , 1 0 ) 
WRITEC6, 22041 

WRITE(6,2205) 

DO  2301  LX  =  1 , NOK 
LWEINC=NOK-LX*l 

WRITE (6,22  02)  ZN(LWEINC) , C T (L WEI NC, KX ) ,K X= 11 , 2 0 ) 

2301  WRI TE (6 , 220  3)  ( T < LWE I NC , KX)  ,K X=2 1, 30 ) 

WRITE (6,2206) 

WRITE(6, 22071- 
DO  2302  LX= 1 , NOK 
LWEINC-NOK-LX*i 

WRITE (6,220  2)  2N( LWE I NC) , < T  (L WEI NC , KX ) , KX  =  31 , 40 ) 

2302  WRITE (6,2203)  <T ( LWE I NC , KX) ,KX=4 1 , 5 0 > 

WRITE (6,2208) 

WRITE(6,2209) 

DO  2303  LX  =  1 , NOK 
LWEINC=NOK-LX»l 

WRITE (6, 2210)  ZN < L WE  INC) , ( T (L HEINC , KX) , KX  =  51 , 5 6) 

2303  WRITE (6,2211)  ( T ( LWE I NC , KX)  ,KX  =  57,62) 

IF(IPRINT)  2401,2401,8000 

8000  WRITE(6,8001) 

00  8101  IQ-1,3 
JO=IQ*IQ-l 

8101  WRITE(b,8021)  ( T  00  A  ( I  Q,  KQ)  ,  KQ=  1, 6)  ,  JQ,  BOD  A(IQ) 

DO  8102  IQ=4 , 6 

JQ= I Q+I Q-7 

8102  WRITE (6,8022)  ( TOOM I Q, KQ ) , KQ= 1, 6) , JQ, BOO A ( IQ) 

WRITE (6,8002) 

DO  8103  IQ= 1 , 3 
JQ=I Q*IQ- 1 

8103  WRITE (6,8021)  ( T OOF ( I Q, KQ ) , KQ  =  1,  6 ) , JQ,  BOOF  ( I Q) 

DO  8104  IQ=4,6 

JQ= I Q*IQ-7 

8104  WRITE<6, 80221  ( T OOF ( I Q, KQ > , KQ=1, 6 ) , JQ , BOOF ( IQ ) 

WRITE (6,80031 

DO  8105  IQ=1,3 
JQ=IQ»IQ-1 
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#105  WRITE (6,8021)  t T008 ( I Q, KQ) , KQ  = 1, 6 ) , JQ, BOOB( I Q) 

00  8106  IQ=4  , 6 
JQ=IQ*IQ-7 

8106  WRITE  (6,8022)  ( TOOB  (  I Q,  KQ)  ,  KQ=1,  6 )  ,  JQ ,  BODBI  IQ) 

WRITE (6,8004) 

DO  8107  IQ=1, 6 

8107  WRITE (6,8023)  BODC( I Q)  ,  ( T OOC ( I Q, KQ) , K Q  =  1 , 6)  , B0D9 { IQ) 

WRITE(6,801i) 

00  8108  IQ  =  1«  3 
JQ=IQ»IQ 

8108  WRITE (6 , 8021)  ( TE VA ( I Q,  KQ) , KQ  =  1, 6) , JQ , BEV A( IQ ) 

00  8109  IQ  =  4» 6 

JQ  =  IQ«-IQ-6 

8109  WRITE (6,8022)  ( TE VA ( I Q, KQ ) , KQ= 1, 6 ) , JQ , BE VA ( IQ) 

WRITE (6,8012) 

DO  8110  IQ=i»  3 
JQ=TQfIQ 

8110  WRITE (6,8021)  ( TE VF ( I Q, KQ ) , KQ= 1, 6) , JQ, BEVE( IQ) 

00  8111  IQ=4, 6 

JQ= I Q+IQ-6 

8111  WRITE(6,8022)  ( TEVE ( I Q, KQ) , KQ=1, 6 ) , JQ , BEVF ( IQ) 

WRITE(6,8013) 

DO  8112  IQ=1,3 
JQ=  I Q  t-IQ 

8112  WRITE  (6,8021)  (TE V0 ( I Q, KQ) , KQ=1, 6) , JQ, BEV B( IQ ) 

00  8113  IQ  =  4 » 6 

JQ=IQ*IQ-6 

8113  WRITE (6,8022)  ( TEVB ( I Q, KQ ) , KQ= 1, 6 ) , JQ, BEVB( IQ ) 

WRITE (6,8014) 

DO  8114  IQ= 1, 6 

8114  WRITE (6,8023)  BE VC ( I Q ) , ( T EVC ( IQ, KQ ) , K Q  =  1 , 6 >, BE VB ( I Q) 

2401  IF  ( IWSTP  .EQ.  1)  CALL  EXCFM 

•ID  M163 
*  D  FH00.30 
*0  FM00.48 
*0  FH00.63 
*0  FM00.82 

•0  FM00.175,FM00.179 

IF(PRNTOP  .EQ.  MIN)  GO  TO  5614 
WRITE (6,2199)  TI TO, H 0 IG l , VK NO TS, EN ( J J ) 

•10  ZMOO 
*0  FMOO . 8 , F H00.9 

970  FORMAT (  F3 . 0 , 5F7. 2 , "5 . 0 , FI 0 .7 ,F5 . 1 , F5 . 0) 

971  FORMAT(12X,F3.0,5F7.2,F5.0,F10.7,F5.i,F5.0) 

972  FORMAT (  4F8.2) 

973  FORMAT ( 12X , 4F 8.2 ) 

•0  FMOO. 12 

7001  FORMAT ( 314) 

•0  FMOO. 20 

2,6HDGAMMA,3X,3HCLZ,5X,3HASP,2X,3HNCW/) 

2215  FORMAT (/15X,4HSWAY,4X ,5 HH E A VE , 3X , 5HPI TCH, 4X , 3HT A W , 4X, 20 HESTIMATED 
2AMPLITU0ES7) 

•0  FMOO. 24, FMOO. 28 

C  HYDROFOIL  VESSEL  WITH  FOIL  POSITION  -  UP=1,  00WN=2 
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C  NUMBER  OF  ZERO  KNOT  CONDITIONS 
C  PRINTOUT  OF  MATRIX  EQUATIONS  -  NO  =  0,  YES  =  t 
REAOJ5*7001)  IFOIL.NZERO.IPRINT 
*D  FM00.59 

C  SPAN,  HYDROFOIL  PIERCING  SURFACE  -  NO=0.  YES=1. 

*  0  FMOO.65 

IF (NZERO  .EQ.  0)  GO  TO  108 
WRITE (6 , 22  IS I 

DATA  CARD  SET  38 


C  ESTIMATED  MOTION  AMPLITUDES  FOR  ZERO  SPEEO  CONDITIONS  GIVEN  BY 

C  VARYING  THE  WAVE  SLOPE  FIRST  FOLLOWED  BY  THE  HEADING  ANGLE.  ORDER 

C  OF  AMPLITUDES  READ  IN  ARE  SWAY  (2) ,  HEAVE(3),  PITCH(5>,  YAW(6) 

C  ( DEG, WORD! 

C 

DO  101  1=1, NZERO 

READ (5, 972)  AMPMOT<2,I) , AMPHOT) 3 , 1) , A MPMOT (5 , 1 ) , AMPMOT  16,1 ) 

101  WRITE (6,973)  AMPMOT ( 2 ,1 ), AMPMOT { 3 , I) , AMPMOT (5 , I ), AMPMOT (6,1) 

108  WRITE (6,2212) 

•O  FM0D.81 

29008(6) ,BEVB(6),T00C(6,6) ,TEVC<6,6» ,B30C(6)  ,BEVC (6) , AMPLI6) 

*1  SP5.59 
NZ=0 
*D  SP5.129 

IF(IFOIL-l)  401,401,400 

400  IF (VKNO TS  .EQ.  0.)  NZ=NZ*1  ”  *  '  "  --  -  -  —  - 

401  DO  999  IWSTP=1,NWSTP 
*0  FM00.128 

C  BEV  MATRICIES.  SUBROUTINE  -ZERO-  IS  USEO  FOR  THE  ZERO  SPEEO  CASE. 
*1  FMOD.132 

IF(VKNOTS)  7200,7102,7104 

7102  AMPMOT (4,NZ) =THM 
DO  7103  JA=2, 6 

7103  AMPL(JA)=AMPMOT( JA,NZ) 

CALL  ZEPO(TOOF,TEVF,BODF,BEVF,GXl,ELL,RHO,NF,Q,TMAS,T,LL,AMPL) 

GO  TO  7106 
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appendix  b 

LISTING  OF  SUBROUTINES  ZERO,  FOIL,  THEO,  and  EXCIT 


The  following  subroutines  are  used  to  calculate  the  motion  coeffi¬ 
cients  and  wave  excitation  forcing  functions  of  the  submerged  hydrofoil 
system  of  a  hull  borne  hydrofoil  craft.  Subroutine  ZERO  is  used  for  the 
special  case  of  zero  craft  speed.  For  speeds  greater  than  zero  but  less 
than  the  critical  "lift-off"  speed,  subroutines  FOIL,  THEO,  and  EXCIT 
are  used. 


SUOROUTINc  7.t  POfroOF,  TEVF.DOOF  ,9E VF  , GX I , E LL , RHO , NF, Q , TMA S  , T , NFREQ , 
ZAMP.W  AM3 L  ,HOG,NOKJ 

COMMON  /Of 01 L/  GA(6,6)  ,EST(5,30, 2)  ,CAL(5,3(J,2)  ,RAT2,RAT6 
DIMENSION  0(10,10) ,AMP(fc)  .CUPTStll) 

DIMENSION  G9(6,6),GC(6,6)  »  OF  (  6 ) , T  13?, 62)  ,  G  A  4  (  6 , 6  ) 

DIMENSION  T00F<6,6)  ,  TEVF I  6 ♦ 6  >  ,  POOF (6)  ,3EVF(6) 

COMPLEX  GO, GC , OF ,  C , VC A , VO ACON 

COMPLEX  C t,C2  ,A3,A4, ARG,HSIM,  MCOS  ,  W  «CCL  ,  WACOM ,  RA.C3  ,WVELL  ,WVELM 
OATA(O^JISlN),N=l,ll)/0.,Q. 04,0.11, 0.31,0. 60, 0,90, 1.13, 1.25, 1,23 
2,1.15,1.07/ 

PlA33-PlA35-PlA22  =  DlA24  =  PlA?6rPlA44=,P2A4  4=P1446=PlA55=PlA66  =  0. 

DO  200  1=1,6 
GE ( I ) =<0.  ,0.1 
OO  201  J  =  1 , 6 

IF(nF°EQ  .EO.  1)  GA(I,J)=0. 

GO  (  I ,  J)  -  (  0  . ,  0  . ) 

201  GC ( I , J)  -  ( 0  .  , 0  . ) 

200  CONTINUE 


r  MULTIPLICATION  FACTORS  FOR  NON-DIM.  ARE 


r  ACCEL  FORCES 

C  VEL.  FORCES 

C  OISPL.  FORCES 

C  SU9SCR I P IS  11,13,31,33,22 

C 

C  INERTIA  MOMENTS 

C  ANGOLA3  VEL.  MOMENTS 

C  ANGULAR  DISPL.  MOMENTS 

C  SUBSCRIPTS  55,44,46,64,66 

C 


< 1 ./MASS) 

C1./MASS)*SQRT(LPP/GRAV) 
C  1  ./MASS)  MLPP/GRAV) 


( 1 ./MASS) /LPP**2 

( l./MASS) *  S  Q  R  T (LPP/GRAV)/LPP**2 
( l./MASS)»(LPP/GRAV)/LPP**2 
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CROSS  INERTIA  < 1  ./MASS ) /LPP 

CROSS  VEL.  I t./MASS) *SQRT  (LPP/GRAV >/LPP 

CROSS  OISPL.  (  J./MASS)  ML  PP/GRAV) /LPP 

SUBSCRIPTS  15, 35,51,53,24,26,42 ,62 

EXCIT.  FORCES/WAVE  AMPL.  L PP/( MASS*GRAV*HAMPL> 
SUBSCRIPTS  1,2,3 

EXCIT . MOM. /WAVE  AMPL.  1 . / t H ASS’GRA V*W AMPL > 
SUBSCRIPTS  4,5,6 

GR A  V=  32 . 2 

RMA  SS  =1 . / TMAS 

ZLOIVG=ELL/GRAV 

ELLSQ  =  ELL  *ELL 

F  A 1 =RMASS 

FA2  =  F Al/ELLSQ 

FA3  =  F  Al/ELL 

F B1 =R  MASS  *SQPT t  ZLOI V  G) 

F  32  =  F  31/ELLSQ 
FR3=F91/ELL 

FOl=ELL/ ( TMAS *GRAV*H AMPL) 

FOZ-F Pl/ELL 

C=CMPLX(0  .,1.) 

PI  =  3. 14159 
RA00EG=57.29573 
XMU= ( 18 1 . -HOG »  7RAOOEG 
SINMUrSIN  tXMU) 

C  0  S  MU  =C  3  S  ( XMU  ) 

Al=°I »RHO 

OMEGA =GX I ’SORT (GRAV/ELL) 

OMEG2=OMEGA*OMEGA 
RK  = OMEGA* OMEGA. GRAV 
awi  =  0 . 5*  A l*OM£GA 


SUMMATIONS  FOP  FOIL  COEFFICIENTS  AND  EXCITATION  TORCES  /  MOMENTS 


00  100  1=1, NF 
CPL=Q  (1,1) 

SPAN=Q( 2 , I) 

CHORD  =Q ( 3 , 1 ) 

X=Q(4,I) 

T=QC5,I) 

Z=Q (5,1) 

OGAMM A=3  (7,1) 

ASP=Q (9,1) 

STRUT  =Q( 13,1) 

MCPL=  CPL 

APE A=SPAN*CHORO 

SPA  N3  =S° A  N*  *  3 

ASPRAT=SP AN/CHORO 

ASPCOR  =  ASPRAT/ (A  SPRAT VASP) 
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CPL=CPL  *  A  SPCOR 
GAHMA=OGAMMA/57.296 
SIN G=SIN( GAMMA) 

SINGSQ=SING*SING 
C OSG= CO S  t  GAMMA) 

COSGSQ=OOSG*COSG 

C  *  *  *  *  *  (FOIL  COEFFICIENTS)  *  *  * 

C  AOOEO  MASS 
C 

IFTNFREQ  .GT.  1)  GO  TO  308 
A2=).  25*AR£A*CHOP.O*CPL 
P1A33=P1A33M  A2*COSGSQ) 

PI A35  =P 1 A  35* ( A2*X*COSGSQ> 

P 1 A  55  =P1 A  55* ( A2*C0SGS0* (CHORD *CHORD/32. +X*X) ) 

P 1 A  22  =P  1 A22  *■  (  A2»SINGS0) 

PlA24=PlA24M  (  A2*  Z  *SINGSQ  )  *  (  A  2  *Y*  SI  NG  *COSG>  ) 

P1A26=P1 A26*( A2*X*SINGSQ) 

PIA44  -PI  A  44*  (A2*  ( ( Z  *S  ING*  Y*C0  SG>*»2>  > 

P2A44=P2A44*(  A  R£  A  *  ARE  A*S  P  AN/4  8  .)  *CPL 
P 1 A46  =Pl A  4b*-  (  (A2*Z*X*SINGSQ>  *  (  A2 * Y*X *  SI NG*COSG>  ) 
P1A66=P1A66* ( A2*SINGSQ* (CHORD • CHORD/ 32 .♦ X* X) ) 

DAMPING 

308  SHAVE  =SW=  0. 

IF(STRUT  .EQ.‘  0.)  GO  TO  38 

SS  =  ( Y*.  5  *  SP AN  *COSG)  *COSG* (Z*. 5*SPAN*SING) *SING 
P=OMEGA»OMEGA*SPAN/GRAV 

li  =  ij  .*=» 

IF ( Ll  .GT.  0)  GO  TO  30 
CW2  =  3  . 

GO  TO  36 

30  IF(Ll-lO)  32,34,34 
32  PI =L1 *1 
®1=P1/13. 

CH2  =  (  10.*  (P-Pl)*  (CUP  TSUI  *2)  -CHPTS(Ll*l>  ) )  *CWPTS(L1*2) 
GO  TO  36 
34  CH2  =  1 .07/P 

36  8HAVE=BW1*CH2*SPAN*AREA*SINGS0*CPL 
38  CONTINUE 

ALF^AeS (93.-OGAMMA) 

MCYCL  =1 

22  IF(AUF-40.)  10,12,12 
10  CD= • 0  467  *  ALF 
GO  TO  14 
12  CO= l . 17 

14  IF(MCYCL-2>  16,17,18 

16  COZ=CO 
MCYCL =2 
ALF  =OGAMM  A 

IF ( ALF  .GT.  90.)  ALF=1S3.-ALF 
GO  TO  22 

17  CO Y=CO 
MCYCL=3 
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S=-Y/Z 

ALF=A8S(57.296*ATAN(S)-0GAMMAI 
GO  TO  22 
18  CO^OLL=0O 
C 

COM  1  =  4. *RH0*CH0R0*CPL/(3.»PIJ 
Ql  =  Y*COSG»-Z*S  ING 

VOA=C*OMEGAM-SING*( AMP ( 2 ) *X* AMP< 6) -Z* AM P I 4» >*COSG*  « A MP< 3 ) ♦ V* AMP ( 4 
2I-X*AHP<5>)> 

VOACON=-VOA 

VOASQR=SQ9T (REAL (VOA*VOACON»  » 

R£LS=SPAN3*REAL(C»0MEGA*AMPt4>  *VOACON>/ ( 12.*VOASQR> 

GB ( 3. 3) =G  0 ( 3 . 3) *COM 1 *COZ * COSG SO* VOA S OR*S P AN 
GFH5,5)=GB(5,5>  «-COMl*COZ*X*X*COSGSQ*VOASQR*SPAN 
IF(STRUI  .ME.  0.)  SW=OWAV£*SINGSQ 
GO  1 2 ,  2)  =GR  ( 2 . 2)  ♦COMl*COY*SINGSQ*VOASQR*SPAN4SW 
I  F  (  S  T  RU  T  .NE.  0.)  SW  =  9WAV£»AB S  (SI MG*01> 

G? ( 4, 2»  =G3(4, 2) - COM  1* COY* SING  *  < VO AS QR*S P AN*Q1+ RELSI 

GO  (4, 4»  =G0(4,4)  «-COMl*COROLL*(  C1M  VO ASQR* SPAN’Ql+RELS)  *RELS*Q1 ♦ « S* A 
2N3*V0 AS3R/12. » J *SW 

GO (4, 6) =G? (4, 6) -COMl*CDY*X*SI MG*C Q1 *  V'„SQR*SPAN *RELS> 

IF(SrR'JT  .NE.  J.)  SW  =  8 WA V  £  * A B S  (X)  *SINGSQ 
GB(6,6)-GB(6,6I fCOM 1 *CO Y* X *X * S INGSO* VO AS QR*SPAN+SW 
C 

c  •  *  *  *  *  (EXCITATION  forcesj  ***** 

50  02=Y* SINNUVX*C0SHU 

Cl  =  -0 .25*PI*RHO* fCH0P0**2» *WA KPL*OMc G2 *C MPLX (COSG , t SINMU*SING>  ) 
C2=CEXP(RK*CMPLX(Zt-Q2) ) 

A3=C1*C2 

A4=RK*CMPLX(SING.(-S I NMU’C jSG ) ) 

A?G=. 5*SPAN»A4 

HSIN=  ,5*(CEXP  (ARGJ-CFXP(-ARG) » 

HCOS=  .5*  (  CEXP  ( ARG) *CEXP (  -  ARG)  ) 

MACCL  =2 • *Cl*C2*HSIN/A4 

W  ACCM  =  ( A  3 / A4 • *(SPAN*HC0S-(2.*HSIN/A4)J 
WA=MAMPL*GV€GA*CMPLX ( ( - S I NMU* S ING ) ,COSGJ 
ARG  =  W  A*C2 

C  3  =  4. *RHO*CHOPO*WA*C2*CSQ?T ( A F G*C ON JG ( ARG) > / < 3 . *PI* A4 ) 

WVSLL=C3  * (SPAN*REAL( A4) *HCOS*2  .* 1 1 .- (REAL (A4J/A4) )*HSTN) 

WVEU1=C3*({2.*Ql-t2.*(l.*R£At(A4>*Ql)/,A4>4.5*SPAN*SPAN*P.EALtA4>M4 

2.*PSAL<44)/(A.*A4»)J*HSTN*-SPAN*(l.»:!EAL(A4>*Ql-{2.*REAL<A4)/A41J*H 

3COSI 

GF (2) =GC ( 2> -( WACCL*WVELL*COYI *  SI  MG* ASPCOR 

GF  (  J)  =Gr  (  3)  ♦ ( WACCL*WVrLL*COZ>  *  COSG*  A  SPOOR 

GF  (  4)  =Gr  (  4)  ♦•  (  (V)ACCU*OlfHACCP)  ♦KVELH  *C  DROLL)  *ASPCOR 

GF (5) =  G(r ( 5  > ♦ C WACCL*WVELL*COZ»  *  x*COS  G* ASPCOR 

GF  (o)  =GF  (5)-  (WACCL*WVElL*COY)  * X*S ING* ASPCOR 

IFIMCPL  .EQ.  1»  GO  ro  10J 
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Y=-Y 

COSG=-COSG 
MCPL=1 
GO  TO  50 
100  CONTINUE 


FOIL  COEFFICIENTS  (NON-OIM.) 


IF(NFREQ 

G  A  ( 3 

,  3)  = 

G  A  (  3 

,5)  = 

G  A  ( 5 

,  3)  = 

G  A  ( 5 

,5)  = 

G  A  ( 2 

,  2)  = 

GA  (2 

,  4)  = 

G  A  ( 2 

.  6)  = 

G  A  (  4 

,  2)  = 

G  A  ( 4 

,  4)  = 

G  A  ( 4 

,  6)  = 

GA  (6 

,  2)  = 

G  A  ( 6 

,4)  = 

G  A  (  6 

,6)  = 

310  GO ( 3 

,3)  = 

GB  ( 5 

,  5)  = 

G  3  ( 2 

»  2)  = 

GB  (4 

.2)  = 

GB  ( 4 

,4>- 

GB  ( 4 

,6)  = 

GB  (6 

,  6)  = 

.GT.  1>  GO  TO  310 


=  FA2M  +A1*P1A55) 


-A1*P1A46> 


00  400  K=l,6 


T (NFREQ 
T ( NFR  fq 
T (NFREQ 
T (NFREQ 
4  j  0  T  (NFREQ 
T (NFREQ 
T  (NFREQ 
T (NFREQ 
T (NFREQ 
T (NFREQ 
T  (NFREQ 
T  (NFREQ 
T (NFREQ 
T (NFREQ 
T  (NFR EQ 
T(NFREQ 
T  (NFREQ 
T  (NFREQ 
T (NFREQ 
T (NFREQ 
T (NFREQ 
T (NFREQ 
T (NFREQ 
T (NFREQ 
T (NFREQ 


,K)=GA(K,K) 

,Kn  J)  =P.EAL  ( GB (K, K ) ) 
,K»2J) =AIMAG(GB(K»K)) 
,K*-30)=REAL(GC(K»tO) 

*  K  *4  1 1 =AIMAG(GC(KrK>» 
, 7  )=GA (3,5) 

,  8 ) =GA  (2* 6) 

,9  )  =  GA (2,4) 

» 1 0 ) =G A (4 ,6) 

,  17)=  PEAL(GB (3,5)  ) 

,  2  7 ) =A INAG (GB (3,5)  ) 
,10)=  REAL (GB (2,6)  ) 
,28) =AIMAG(G9 (2,6)  ) 
,19)=  REAL (GB (2,4)  ) 
,29)=AIMAG(GB (2,4) > 
,20)=  REAL(GB (4,6)  ) 

,3 ))=AIMAG(GB (4,6)  ) 
,37)=  REAL(GC (3,5)  ) 

, 4 7 ) = A IMAG (GC (3,5)  ) 
,38) =  REAL (GC (2,6)) 
,48)=AIMAG(GC (2,6)) 
,39)=  REAL(GC  (2,4)  ) 
,49) = A IMAG(GC (2,4) ) 
,40)=  R£AL(GC (4,6)  ) 
,53)=AIMAG(GC  (4,6)) 
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MULTIPLICATION  OF  ACCELERATION  TERMS  BY  -GXI*GXI 
MULTIPLICATION  OF  VELOCITY  TERMS  BY  *GXI 


A=-GXI*GXI 

B=GXI 

00  232  1=1,6 
00  203  J=l*6 
G A A (I ,J) =A*GA (I,  J) 

203  G8(I, J)=B*G3(I,J> 

202  CONTINUE 

FOIL  COMPONENTS  FOR  MATRICIES  -TOO  AND  TEV- 
T00(R0N,C0L)  ,  T  EV  ( R  CM, COL) 

PRINTOUT  OF  MATRICIES  IS  FOR  THE  FIRST  WE(ND)  PRINTED 

00  205  1=1,3 
00  236  J= 1 , 3 
I£V  =  I *1 
JEV=JvJ 
100=1 EV-1 
JOD  =  J  EV- 1 

TODFl I, J) =GAA (IOO,JOO>-AIMAG{Ge<IOO,JOO) ) +RE AL ( GC (IOD , JOO) ) 
T  OOF  (  I,J+3)=*R£AL(GB(I00»J00)  )  *A  I  MAG  ( GC  (  1 00 ,  JOO )  ) 

TOOF( 1+3 ,  J  +  3 ) =TOOF(I,J) 

T00F(I+3, J)=-TOOF(I, J^3) 

TEVF( I , J ) =GAA (IEVtJEV) - A  I  MAG ( Ge ( I EV , JEV ) ) ♦REAL (GCIIEV, JEV)> 
TEVF( I, J+3) =f REAL<G8(I£V, JEV) ) ♦AI MAG ( GC ( IE V , JE V > ) 

TEVF(  I  ♦  3  ,  Jv3)“=TEVF  Cl  ,J) 

TEVF( 1*3, J) =-TEVF(I, JM3) 

206  CONTINUE 
2C5  CONTINUE 


EXCITATION  FORCES  AND  MOMENTS  (NON-DIM.) 


GF ( 2) =F  0 1 *GF (2) 

GF(3)=F01*GF<3) 

GF(4)=FD2*GF(4> 

GF(5) =F02*GF(5) 

GF (6) =F02*GF (6) 

DO  412  <=1,6 

T(NFR£Q,K>5))=  REAL(GF(K)) 

432  T (NFR=Q,K*56> =AIMAG(GF (K> ) 

FOIL  COMPONENTS  FOR  MATRICIES  -BOO  AND  BEV- 

BOOF ( 1) =3  .0 
800F(2)=»RErtL(GF(3)) 

BODF ( 3) =»REAL (GF(5)) 

BOOF ( 4) =0.0 

ROOF ( 5) =  - AIM A  G ( GF ( 3)  ) 

800 F( 6) =- AIMAGCGF(5> ) 

3EVF(  1) =fREAL (GF (2) ) 

8EVF ( 2) =  fREAL (GF (4) ) 

8EVF( J) =>REAl  (GF (6 ) ) 

8EVF  ( 4)  =  -AMAG(GF(2)) 

BEVF(5)=-AIMAG(GF(4) ) 

9FVF ( 6) =-AIMAG(GF(6) ) 

RETURN 
999  ENO 
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SUBROUTINE  FOIL(TOOF,TEVF,BODF,9EVF,VKNOTS,WAMPL,HOG1,GXI,ELL,RHO, 
2NF  ,  Q,  TMAS  , T , NFREQ, W AVL ) 

COMMON  /QFOIL/  GA16.6) 

DIMENSION  T00F(6,6) ,TEVF<6,6>  , OOOF 1 6) ,3E VF C 6) , Q < 1 3 , 10 ) 

DIMENSION  GO ( 6, 6 ) ,GC  (6,6)  ,GF | 6  ) ,T (32 . 62) ,GA A IS , 6 ) 

COMPLEX  CK,B2,C2,P1B33,P1C33,  P2B3 5 , P1C35 , P2C35  , P1B22, P1B24,P1C24 ,  P 
12826,P1C26,P1B44,P2R44,P1C44,P2946,P1C46,P1B53,P1C53,P3B55,P1C55 ,P 
!  22C55,P1B62,P1B64,P1C64,P2B66, P1C6S 

COMPLEX  C ,EXL,EXM,AA ,ARG,HSIN  ,  HCOS, VI , V2 , W1 , W2 , PL1, PF 2, PF3, PL2,  3F 4 
2,PM1,PF5,PF6,88,W3 

,  COMPLEX  GB.GC.GF 

%  P1A33=P1A35=P1B35=P1A22=P1A24=P1A26=P1B26=P1A44=P2A44=P1A46=P1B46= 

1P1A55=P1B55=P2B55=P1A66=P1Q66=P3366=0. 
P1833=P1C33=P2835=P1C35=P2C35=P1B22=P1B24=P1C24=P2B26=P1C26=P1844= 
Is  2P2844=P1C44=P2B46=P1C46=P1B53=P3B55=P1C55=P2C55=P1B62=P1B64=P1C64= 

.  3P2B66=P1C  66= (  G  .  ,  0  • ) 

PF2=PF3=PF4=PF5=PF6=<0.,0.) 

00  230  1=1,6 

i  GF(I)=I0.  ,0.1 

DO  231  J  =  1 , 6 

IF(NFR£3  .EQ.  1)  GA(I,J)=Q. 
jp  G3(I, J)=<0.,0.) 

231  GC(I,J)=(3.,3.I 

233  CONTINUE 

MULTIPLICATION  FACTORS  FOR  NON-OIM.  ARE 
ACCEL  FORCES  II, /MASS) 

VEL.  FORCES  I  1  ./MASS) *SQRT  (LPP/GRAV) 

OISPL.  FORCES  I l./MASS)*(LPP/GRAV) 

SUBSCRIPTS  11,13,31,33,22 

INERTIA  MOMENTS  (  1  ./MASS) /LPP**2 

ANGULAR  VEL.  MOMENTS  I  1  ,/M AS S) *SQR T ( LPP/GRAV ) /LPP* *2 

ANGULAR  OISPL.  MOMENTS  ( 1  ./H ASS >* I LPP/GRAV ) /LPP* * 2 
SUBSCRIPTS  55,44,46,64,66 

• 

CROSS  INERTIA  I  1  ./M ASS ) /LPP 

CROSS  VEL,  (  1  ./MASS) *SQRTCLPP/GRAV)/LPP 

CROSS  OISPL.  I  l./MASS)* (LPP/GRAV ) /LPP 

SUBSCRIPTS  15,35,51 ,53,24,26,42 ,62 

EXCIT.  FORCES/WAVE  AMPL.  L FP/ ( MA SS*GR AV*WAMPL ) 

SUBSCRIPTS  1,2,3 


EXCIT. MOM. /WAVE  AMPL.  1 ./( M ASS  * GRAV*W AMPL ) 

SUBSCRIPTS  4,5,6 


o  o  o  o 


GRA V=32 • 2 
RMASS=1 . / TMAS 
ZLDIVG=ELL/GRAV 
ELLSQ=ELL*ELL 
F  A1=RMASS 

F01=RMASS*SQRT(ZL0IVG) 

FC1=RMASS*ZL0IVG 

FA2=FA1/£LLSQ 

FB2-FR1/ELLSQ 

FC2=FC1/ELLSQ 

FA3=F Al/ELL 

FB3=F01/ELL 

FC3=FC1/ELL 

F01=ELL/ ( TMAS*GRAV*HAMPL) 
F02=FD1/£LL 

PI=3. 14159 

XMU=(  183.  3-HOGl)/57. 2  95  77  95 
U=VKN0TS*1.689 
S INMU  =S IN  (XMU) 

COSMU=COS  fXMU) 

OMEGA =S3RT (2. *PI*GRAV/WAVL) 

OMEGA £  =  3XI*SQPT  (GRAV/ELL) 

Al=PI*RHO 

A3=A1*U 

81=9. 5»RHO*U 

C 1=  01 *U 

BB=CMPLX(0.,1.) 


SUMMATIONS’  FOR  FOIL  COEFFICIENTS  AND  EXCITATION  FORCES  /  MOMENTS 


00  13  0  I  =  1 >  NF 
CPL=0<1.I1 
SPAN=Q(2, I) 

CHORO  =Q ( 3  » I) 

S  =  0(4,U 
r  =  Q(5,Il 
Z  =  Q  (6 , 1 ) 

OGAMMA=a(7,I) 

CLZ=Q(8,I) 

ASP=Q  (9,1) 

NCPL=CPL 

ARE A=SP4N»CHORO 

ASPRAT=SPAN/CHORO 

ASPCOR=ASPRAT/( A  SPRAT* ASP) 

CPL=CPL*ASPCOR 

G A MM A =0 GAMMA/ 5 7. 2957795 

SING=  SIN ( GAMMA) 

SINGS 0=SING*S ING 
COSG=COS (GAMMA) 
COSGSQ=COSG*COSG 
SI  N*1UG  =  S  IHMU*  SING 
CL ALP  H=  2 . »PI 
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XK1  =  0  .5*omegae*choro/u 
XK2=(0MEGA*0H£GA)/GRAV 

XK3=0 .5* (0MEGAE*0MEGAE»  *C HORO *COSMU/GRA V 
CALL  THEO(XKl.CK) 

A2=  0 . 25*AREA*CHOP.O*CPL 

82=AR£A»CLALPH*CK*CPL 

C2=CLZ*ARFA*CK*CPL 

C*  *********** 

IFCNFPE3  .GT.  1)  GO  TO  308 
PlA33=PlA33+(  A2*C0SGSQ) 

P1A35=P1 A35+( A2*S*C0SGSQ> 

Pi A55  =P1 A  55  + (  A2*C0SGSQ* (CHORD  ‘CHORD/  32.+S*S) ) 

PlA22=PlA22+( A2*SINGS0) 

P1A24=P1A24+ (  (A2*Z*SIHGSQ)  ♦  < A 2*Y *  SI NG *COSG>  1 
PI 4  25  =P 1 A  26+ ( A2*S*Sr  NGSQ) 
P1A44=P1A44+(AR£A*AR£A*SP4N/48.)*CPL 
P2 A  44 -P 2 A  44 K A2* ( (Z* SING  +  Y*COSG> **2) ) 

PlA46=Pl A  46+ ( (A2*Z*S*SINGSQ) ♦ ( A2+ Y* S*SI NG*COSG > ) 

PlAo6=PlA  66  +  ( A2  +  SINGSO* (CHORO  +CH0R0/32. +  S*S>  > 

3C8  Pl333=PlB33+(«2*COSGSQ> 

P1B35=P1B35+ ( A2*COSGSQ) 

P2835=P2B35+ ( B2 *COSG SO  * ( S  +  (CHORD/ 4. ) ) ) 

P1353-P1B53+(B2» (S- ( CHO RO/4 . > )*COSGSQ> 

P 1355=? IB  55+ ( A2*S*COSGSQ) 

P2355=P2855+( (CHORD* *31 *SF AN* COSG SO/ 16 . )  *CPL 
P2B55=P)B55+ (02* (S+  (CHOP.O/4.  )  )  *  (S -  ( C HORD/4  .  )  )  +COSGSQ) 
P1822=P1B22+(B2*SINGSQI 

P 1 024  =P1 B24* ( (32*Z*SINGSQ )+(82*Y*SING*C0SG) ) 

P1825=PlB26+( A2*SINGSQ) 

P?926=P2H26+ ( B2* (S+ (CHO R 0/4. ) ) *SINGSQ> 
P1944=P1044+(P2*SPAN»SPAN/12. I 
P2944=P2044+ (H2* ( (Z*SING+ Y*COSG) **2)  ) 

P1B  46=P 1 B  46  + (  (A2*Z*SINGSQ) ♦  ( A  2 *Y* S I NG*COSG> ) 

P2B46=P2B46+  (B2* (S  + (CHOPO/4. >  ) * ( t Z» S I NGSQ»  + C Y*S X NG*CO SG> ) ) 
P1862=P1062+(P2* (S- (C H0R0/4 . >  )*SXNGSQ) 

P1364=P1R64+ (02* (S- ( CHORD /4. > ) * M Z* SI NGSQ) + 1 Y* S I NG*COSG> ) ) 
P1366=PlB66+( A2*S*SINGSQ) 

P2366=P2066+  102* (S+ ( CH0R0/4 . )  t  *( S- ( CHORD/4 . )  )  *SINGSQ> 
P3366=P3B66+ ( (CHORD  *  * 3 ) *S PAN*  SING  SO/ 16. ) *CPL 
P1C33=P1C33+(C2*C0SG) 

P1C35=P1C35+(02*C0SGSQ) 

P2C35=P2C35+(C2* ( S- ( CHORO/4 . ) ) *C0SG) 

P1C53=P2C35 

P 1C55 =Pl C 55 +(02*(S-( CHORO/4.) ) *C0SGSQ) 

P2C55=P2C55+(C2*  (S  +  (CHOP.D/4.  »  >  *(  S - ( C H0R0/4 . >  )  *C0SG> 
P1C24=PIC24+(C2*Y*SING> 

P1C26=P1C26+(B2*SINGSQ) 

P1C64 -PIC  64+ (C2»  <S- (CHORD /4. )) *Y*  SING) 

PlC66=PlC66  +  (B2*(S-(CH0R0/4.)  )*SINGSQ)  • 

P 1C  44 -P 1C44+ ( C2*Y* ( (Y*COSG)  +  ( Z*SING) ) ) 

P1C46=P1C46+ ( ( 02  *Z*S INGSQ > ♦ ( B 2 *Y *  SI NG*COSG) ) 
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iri 


o  o  o  o 


c  * 


************ 

CALL  EXCIT(XK3,XK1,CK,EXL,EXM> 

25  C  =  CMPl_X  (COSG,  SINMUG) 

XRE  AL=X<2*SING 
XIMAG=X<2*SIMMU*COSG 
4A=CMPLX(  XREAL*-XIMAG) 

APG=0 . 5*  A  A*SP AN 

HSI N=  3 • 5  *  (CEXP(AP.G)  -CEXP(-AfiG>> 

HCOS=  3.5* (CEXP1ARG) *CEXP (-ARG »> 

V 1  =  (2  ./A  A  )  *HS  IN 

V2=(1./(AA*AA))  * ( AA*SPAN*HC0S-2.*HSIN) 

XREAL  =X  <2  *Z 

XIM AG=-XK2*Y*SINMU 

AA=CMPLX (XREAL.XIMAG) 

AA=WAMPL*OMEGA*CEXP( AA)*C 
Wl= AA*V1*88 
W2=AA*V2*9B 
W  3  =  A A  *08 

IF(NCPL-l)  51,51,50 
53  Wl=2. *H1 
W  2=  2. *W 2 
W3  =  2. *W3 

51  X  I M  AG  =- XK  2*S*  COSMU 

AA=CEXP(CMPLX (0. , XIM AG) ) 

C=CHORO»EXL*A  A 

PL1=C‘AIMAG(H1) »BB 

PF  3=P  E3* (PL1*CQSG*  ASPCOR) 

PL1=C*REAL(W1» 

PL2=CH0R0*W2*EXL 

PM 1=0 •25*CHORO*CHORD*EXM*AA 

PF5  =PK5 - ( I <-S*PLl)-PMl*88*AIMAG(W3)  I  *COSG> *ASPCOR 
HQ=MDG1 

I F ( HO  .3T.  183. »  HO=HP-183. 

I F ( HO  .GF.  172.)  GO  TO  100 
IF  (HO  .'-T.  8.)  GO  TO  100 
PF2=PF2- (PLl*SING*ASPCOR> 

PF4  =  PF4  *•  (PL2*PLi*(Y*COSG*Z*SING))  *ASPCGR 
PF6=PF6«-  (  ((S*PL1)*PM1*REAL(W3))*SIN&»*ASPC0R 
100  CONTINUE 


FOIL  COEFFICIENTS  (HON-OIM.) 


IF ( NFRE3  .GT.  1)  GO  TO  310 
GA( 3, 3) =FAt*( *AI*P1A33) 

GA( y, 5) =FA3*(-A1*P1A35) 

G  A  (  5,  3)  =  G  A  ( 3 , 5) 

GA15,5)=FA2* (*AI*P1A55> 

GA(2,2)  =FA1M  *A1*P1A22> 

GA(2,4) =FA3*(-A1*P1A24) 

G  A ( 2 , 6) =FA3*( +A1*P1A26) 

G  A  (  4,  2J=GA(2,4J 

GA(4, 4) =F A2* ( v A1 * ( PI A44 FP 2 A44  > ) 
G  A ( 4, SI =FA2*(-A1*P1A46) 
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GA(6,2)=GA(2,6>. 

GA (6, 4) =G A (4.6) 

GA (6,6) =FA2*( «-Al*PlA66) 

310  GB(3,  3)  =FB1*(  4-fU*Pin33> 

GB(3,5)=FB3M-A3»P1B35-01*P2B35> 

GB(5, 3) =F33*(-01*P1B53) 

GB ( 5 , 5) =FB2* ( ♦  A 3 »P1 85 5* A 3 * P28  5 5«-B  1*P 3B55  ) 
G9(2,  2)  =FB1*U01»P1B22> 

GB(2, 4) =Fn3*(-01*PlB24) 

GB(2,6) =  FB3  * ( ♦A3*P1B26»B1*P?8  26) 

GB (4, 2) =GS(2, 4) 

GB(4»4)=F82*(»R1MP1E-,4*P2B44)) 
G8(4,6)=F<?2M-A3  »P 1 B 46- B 1 • P2 3 4  6) 

G8(6, 2) =FB3»<  H)1*P1362) 

GB (6, 4) =FB2*(-B1*P1064) 

G3(6,6J=F82*(  *•  A  3  *P  1  B66*  B1  »P2B  E6*  A  3*P3  B66  ) 
GC(3,3)  -FC1M  ♦C1*P1C33) 

GC( 3,5) =FC3*( -C1*(P1C35»H2C35 )  ) 

GC (5, 3) =FC3» ( -C1*P1C53) 

GC(5,5)=FC2M fCl* (P1C55*P2C55  )  ) 

GC  (2,  4)  =FC3M  -C1*P1C24) 
GC(2,6)=FC3M»C1*P1C26) 

GC  (4,  4)  =  FC2M  *-Cl*PlC44) 

GC( 4, El =F  C2* (~C1*P1C46) 

GC(6. 4) =FC2*(-C1*P1C64> 

GC(6, 6) =FC2*( tCl *P1C66) 

00  4 3  G  <=1,6 
T (NFREQ, K)=GA  (K,K) 

UNFREQ.KMO) =REAL (GB(K,K) ) 

T(NFREQ,K*2])  =  AIMAG(GB(K,K)I 
T (NFREQ, K  +  33) =REAL (GC (K,K>  > 

433  T  ( NFR  EQ ,  K  »4  3 )  -  A  I  MAG ( GC ( K, K) ) 

T (NFREQ, 7)=GA  (3,5) 

T (NFREQ, « ) =GA  (2,6) 

I ( NFR  EQ, 9 )  =  G A (2,4) 

T (NFRtQ. 1 3) =GA  (4,6) 

T (NFREQ ,17)=  RE A L < GB ( 3 , 5 ) ) 

T( NFREQ, 2 7) =AIMAG(GP  (3,5) ) 

T(NFREQ,13)=  REAL (GO  (2, 6 ) ) 

T (NFREQ,28)=AIMAG(G8 (2,6) ) 

T  (NFREQ»19)=  PEALJGB  (2,4)  ) 

T (NFRSQ,29)=AINAG(GB(2»4)) 

T  (NFREQ, 2  3)  =  REAUGB  (4,6)  ) 

T (NFREQ, 3  3) =AIMAG(GS  (4,6) ) 

T (NFREQ, 37)=  RE AL ( GC ( 3 , 5 ) ) 

T (NFRtQ, 47)=AIMAG(GC (3,5) ) 

T (NFREO, 3  8)  =  PE  A L ( GC (2 , 6 ) ) 
T(NFREQ,4fl)=AIPAG(GC  (2,6) ) 

T (NFREQ, 39)=  REALJGC  (2,4)  ) 
MNFR,EQ,49)=AIMAG(GC(2,4>  > 

T (NFREQ ,43)=  PE AL (GC (4,6 ) ) 

T (NFREQ, 53)=AIMAG(GC  (4,6) ) 
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ooo  noon  o  o  o  o  o  o  o 


MULTIPLICATION  OF  ACCELERATION  TERMS  BY  -GXI’GXI 
MULTIPLICATION  OF  VELOCITY  TERMS  BY  +  GXI 


A  =  -GX 1*3  X  I 
8=GXI 

DO  202  1=1,6 
DO  203  J  =  1 » 6 
GAA (I ,J)=A*GA  (I,J) 

203  G9(I, J)=B*G3(I,J> 

202  CONTINUE 

FOIL  COMPONENTS  FOR  MATRICIES  -TOO  ANO  TEV- 

00  205  1=1,3 
OO  206  J= 1, 3 
IEV=I +  1 
JE V  =J  +J 
100= I FV- 1 
JOD  = J  EV -  1 

TOQF( I, J) =G A A ( IOD , JOO) -A  I  MAG t  G BC 1 00 , JOO)  >  f REAL ( GC (XOO , JOO) I 
T00F(I,J+3J=+REAL(GB  tlOO, JOO)  )  +AIMAG (GC ( IOO, JOO) > 

T00F( 1+3 . J  +  3) =TOOF ( I , J ) 

T OOF ( 1  +  3  «  J) =- TOOF ( I  «  J  +  3) 

TEVF(  I,  J)  =GA  A  (IEV.JEV)  -  A  I  MAG  (  G9t  I E  V ,  JEV)  )  4-REAL  ( GC  ( I  EV  ,  JEV)  1 
TE'/FI  I,  J  +  3)  =+REAL(G9  <IEV,  JEV)  )  +AIM AG (GC ( IE V , JE V >  ) 

TEVF( 1+3, J  +  3) =T£VF ( I , J) 

TEVF( 1+3, J)=-TEVF(I, J+3) 

2C6  CONTINUE 
205  CONTINUE 


EXCITATION  FORCES  ANO  MOMENTS  (NON-DIM.) 


GF(2)=FD1*(+A3*PF2) 

GF(3)=F01M  +  A3*PF3) 

GF  ( >♦ )  =F0  2  M  + A  3*PF4) 

GF(5)=F02M+A3*PF5) 

GF(S)=FD2M+A3*PF6) 

00  402  <=1,6 

T  (NF.REQ*  K  +  50)  =  RE  AL  (GF  (K  )  ) 

402  T (NFR£Q,K*56) =AIMAG(GF(K> J 

FOIL  COMPONENTS  FOR  MATRICIES  -BOD  ANO  BEV- 

B00F( 1) =0  .0 

900<M  2) =  +  REAL (GF(3) ) 

BOOF ( 3) =+REAL (GF(5)) 

900F ( 4) =  0  . 0 

BCOF ( 5) =-AIMAG(GF(3) ) 

BOOF ( 6) =  - AI M A  G ( GF (5) ) 

BEVF ( 1) = ♦ RE  AL (GF ( 2) ) 

3EVF(2)=+REAL(GF(4)) 

9EVF( 3)  =  +REAL  (GF (6) ) 

BEVF (4) =- AIMAG(GF(2)» 

BEVF(5) =-AIMAG(GF(4) I 
9E VF ( 6) =-AIMAG(GF(6) ) 

RETURN 
999  ENO 
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DTNSROC  ISSUES  THREE  TYPES  OP  REPORTS 

DTNSROC  REPORTS.  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY.  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


